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ABSTRACT 



We present high precision photometry of Kepler-41, a giant planet in a 1.86 
day orbit around a G6V star that was recently confirmed through radial veloc- 
ity measurements. We have developed a new method to confirm giant planets 
solely from the photometric light curve, and we apply this method herein to 
Kepler-41 to establish the validity of this technique. We generate a full phase 
photometric model by including the primary and secondary transits, ellipsoidal 
variations, Doppler beaming and reflected/emitted light from the planet. Third 
light contamination scenarios that can mimic a planetary transit signal are sim- 
ulated by injecting a full range of dilution values into the model, and we re-fit 
each diluted light curve model to the light curve. The resulting constraints on the 
maximum occultation depth and stellar density combined with stellar evolution 
models rules out stellar blends and provides a measurement of the planet's mass, 
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size, and temperature. We expect about two dozen Kepler giant planets can be 
confirmed via this method. 



Subject headings: tbd 



Introduction 



NASA's Kepler satellite has been photometrically monitoring more than 150,000 main- 
sequence stars since its launch in 2009. The primary goal of the mission is to determine 
the frequency of Earth-size planets in the habitable zones of Sun-like stars, and in this 
uest 2326 planetary candidates have been identified with the first 16 months of flight data 
Batalha et al.ll2012l ). Of these, 203 are giant planets with radii (R p ) between 6-15 times the 
radius of Earth (-R©). With Kepler's high photometric precision, both the primary transits 
and secondary transits (occultations) can be measured for many of these giant planets. 

Occultation measurements allow us to better characterize planets by providing con- 
straints on the size and orbital parameters of a companion that can produce the shape of 
the light curve. In addition, depending on the wavelength at which they are gathered, they 
can provide information to estimate the effective temperature and reflectivity of the planet. 
The use of phase curves (the variations in the light curves of a star+pla net system as the 
planet orbits the star) as a means to detect exoplanets was presented by IJenkins &: Doyle 
( 120031 ) . Their analysis predicted hundreds of close-in giant planets with periods up to 7 days 
could be detected by Kepler via their reflected light signatures. More recent studies have 
shown that both transiting and non-transiti ng planets can be detected by measuring the 



variations in light induced by the companion (IMazeh fc Faiglei 



Shporer et al.ll201ll ; iMislis fc Hodgkinll2012l ; IMazeh et al.ll2012l ). Herein, we present a new 



20ldlFaigler & Mazehll2011 



method to confirm giant planets based solely on the analysis of light curves by modeling 
both transits and occultations and eliminating other potential non-planetary sources that 
could produce the shape of the light curve. 

Nearby stars that are captured within the target star aperture can dilute the transit 
signal, resulting in an underestimate of the transit depth and therefore the size of the tar- 
get star's companion. These third-light or 'blend' scenarios can include a background or 
foreground eclipsing bina ry star system or a physically b ound stellar companion in a hierar- 
chical triple star system ( jSeager fc Mallen-Ornelasll2003l ). each of which has the potential to 
produce a transit-like signature. 



Confirmation of Kepler candidate planets typically requires additional ground-based 
follow-up observations due to the prevalence of astrophysical false positives. These tech- 
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niques, which include spectroscopy, speckle and adaptive optics imaging, precise Doppler 
measurements and combinations thereof, are used to help eliminate blend scenarios in order 
to confirm that a Kepler planetary candidate is indeed a planet. The number of planets 
that can be confirmed in this manner, however, is limited by the availability of telescope 
time. In this article, we present a new method of confirming giant planets without the 
need for follow-up observations. We generate a full phase photometric model light curve 
that includes the primary transits, occultations, ellipsoidal variations, Doppler beaming and 
reflected/emitted light from the planet. We then inject a full range of dilution values into 
the model to simulate third-light contamination, then re-fit the diluted light curve models 
to the photometry. Comparison of the fitted parameters with stellar evolution models can 
eliminate systems that are inconsistent with a stellar blend. 

To demonstrate this confirmation method, we analyzed the photometry of a star in 
the Kepler Field of View which shows a signature of a transiting Jupiter-sized planet in 
a 1.86 day orbit. Kep ler-41b was recently confirmed through radial velocity measurements 



(ISanterne et al.ll2011al ). thereby providing a good test case for our method. In Section 2 we 
present the photometry of Kepler-41 and discuss our method to correct for systematics and 
stellar variability. Our full phase photometric model and best fit parameters are presented 
in Section 3, and the confirmation technique and results are discussed in Section 4. Albedo 
estimates for Kepler-41b and other hot Jupiters are discussed in Section 5, and Section 6 
provides a summary. 



2. Kepler Photometry 



Kepler-41 (RA=19 /l 38 m 03M8 and Dec=45°58'53.9"), also iden tified by Kepler id enti- 
fication number (KID) 9410930 in the Kepler Input Catalog (KIC) farown et al-lboilh and 
by Kepler Object of Interest number KOI-196, is a G6V star with an apparent magnitude in 
the Kepler bandpass of K p = 14.465. Kepler-41b was id entified as a Jupiter- sized candidate 
companion to this star in Kepler's second data release ( IBorucki et al.ll201ll). The SOPHIE 
spectrograph obtained radial velocity measurements of this object (ISanterne et al.ll2011al ). 
which led to an estimated mass of 0.55 ± 0.09 Mj. 

The Kepler observations of Kepler-41 described in this article were acquired between 
2009 May 13 and 2011 March 5 and include quarters Ql through Q8. Data were sampled 
nearly continuously during each quarter at 29.42 minute long cadence (LC) intervals (where 
each LC includes 270 summed 6.5 second exposures). Although short cadence data (which 
are sampled more frequently at 58.85 s intervals) are more sensitive to the ingress/egress of a 
transit and can provide better constraints on the mean stellar density, they are not necessary 
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to measure the phase curve and also are not away s available. The raw pixels collected for 
these stars were calibrated (jOuintana et al.l 120101 ) and aperture photometry, background 
removal and cosmic ray corrections were performed with the Photometric Analysis (PA) 
software maintained by the Kepler Science Operations Center to produce the light curves 
( jTwicken et al.ll2010l ). Note the light curve data for this target are publicly available at the 
Mikulski Archive for Space Telescopes (MAST0). 

The effects of instrumental signals in the flux time series were mitigated by fi t ting; an d 
subtracting cotrending basis vector^ (their use is documented in Barclay et al. ( 2012a )) 
from the light curve using the PyKE softwar^E We used the first four cotrending basis 
vectors and fit them to the data using a linear least-squares approach. The light curve 
was stitched together by normalizing the flux by the median value per quarter. Outliers 
and additional remaining signatures due to instr umental arti f acts ( such as those due to 
'sudden pixel sensitivity dropouts', as described in ISmith et al.l ( 120121 ) ) were then identified 
and removed. In total, 215 measurements were discarded yielding a total of 29,680 LC 
measurements. 



We next applied a Fourier decomposition algorithm to separate out star spot-induced 
variability from the light curve. Assuming a coherent signal due to the planet and a constant 
orbital period, sinusoidal components to the light curve were iteratively fit and removed 
to filter out all frequencies that were not affiliated with the planet orbital period and its 
associated harmonics. Specifically, if the amplitude of the peaks in the Fourier Transform 
were >3.6 times the standard deviation (corresponding to 3a), we removed those frequencies 
from the light curve. The PA-corrected, cotrended, and Fourier-filtered light curves for 
Kepler-41 are shown in Figure 1. The next section describes our light curve modeling to 
compute the best-fit planet parameters. 



3. Model Fitting 



Our full ph ase photometric model uses a circular orbit (eccentricity = 0) and we adopt 
the formalism of iMandel fc Agoll (120021 ) to compute the light curve model. Figure 2 shows the 
phase-folded and binned photometry for the primary transit (lower curve) centered at orbital 



1 http: / / archive.stsci.edu/kepler 

2 The cotrending basis vectors are available from http://archive.stsci.edu/kepler/cbv.html 

3 PyKE is a suite of Kepler community 
http : / /keplergo . arc . nasa . gov/PyKE . shtml. 



software 



available 



from 
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phase = with the best-fit model shown in red. Here, <fi = 2n(t — T^/P, where T$ is the 
epoch (the time of first mid-transit) and P is the orbital period of the planet. The occultation 
(top curve) near <fi = 0.5 has been magnified (see top and right axes) and the best-fit model 
is shown in green. Our model of the transit for Kepler-41 includes nonlinear limb darkening 
with four coefficients that we compute by performing a trilinea r inter polation over T e g t , logg, 
and Z = [Fe/H] using tables provided bvlClaret fc Bloemenl (120111) . The stellar properties 
(T efft , \ogg, and Z) were adopted from ISanterne et al.l ( 12011al ) (shown in Table 1), and the 
limb darkening coefficients we use are [cl, c2, c3, c4] = [0.6325, -0.3162, 0.8937, -0.4240]. 
The occultation is modeled in the same manner but we assume the companion is a uniform 
disk (we neglect limb darkening) due to the relatively short ingress and egress times. Our 
model of the phase-dependent light curve takes into account the photometric variability that 
is induced by the companion, which includes ellipsoidal variations (F e n), Doppler beaming 
(F dop ), and light contributed by the planet (-F ref , which includes both reflected star light and 
thermal emission). The relative flux contributions from each of these time-dependent effects 
are distinct and can be decomposed as 

F* + F e ii(t) + F dop (t) + F ie{ (t) 



F(t) 



tot 



where F* is the illumination measured at phase </> - 
from the companion, and F tot = + F ref (0 = 0). 



1 



0.5 when the star is blocking the light 



Ellipsoidal variations in the light curve are caused by changes in the observable surface 
area of the star due to tidal distortions induced by the co mpanion (IPfahl et al.ll2008l) . They 
have previously been detecte d in eclipsing binary s t ars (Wilson & Sofia Il976h~ and more 



recently in exoplanet systems (IMazeh fc Faigler 2010 ; Faigler fc Mazeh 2011 ; Shporer et al. 



201ll ; IWelsh et al.ll2010al : Mazeh et al.ll2012l : Mislis fc Hodgkinll2012l ; Mislis et al.ll2012f ). The 



amplitude of the ellipsoidal var iations is roughly equal to the ratio of the tidal acceleration 
to the stellar surface gravity (IPfahl et al.l 120081 ) assuming tidal equilibrium, and can be 
approximated by 



A 



ell 



OL e ii—— I 

ell M | 



a 



sin 2 i 



(2) 



where M p and M* are the masses of the planet and star, respectively, i?* is the stellar radius, 
a is the semimajor axis, i is the inclination of the orbit relative to the line-of-sight. The 
parameter a e u is defined as 

0.15(15 + «)(! + g) 



OCdl 



(3 



u 



(3) 



wher e u and g are the limb darkening and gravity darkening coefficients, respectively ((Morris 
19851 ). We compute u = 0.6288 and g = 0.4021 by linearly interpolating over T ef[it , logg, and 
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Z using tables provided by lClaret fc Bloemerj (1201 lh . For a circular orbit, the contribution 
of flux from ellipsoidal variations, which oscillate on timescales of half the orbital period (see 
Figure 3), is 



-A e ll COS(27T0 X 2) 



(4) 



tot 



Doppler beaming is an apparent increase/decrease in stellar flux due to Doppler shifts 
in the stellar spectrum that are caused by the reflex star motion around the center of mass 
due to the companion. These signals have on ly recently been measured in transit light curves 
(jMazeh &: Faiglerll2010t IShporer et al.ll20! lh and oscillate with the orbital period. Note that 
we did not detect variations from Doppler beaming in the Kepler-41 light curve, but we 
include a description here because it may be applicable to other planet candidates. The 
amplitude of this signal is 



.4, 



dop 



-a dop AK 



(5) 



where c is the speed of light and ad op is a Doppler boosting factor that depends on the 
wavelength of observation and o n the stellar spectrum . We compute ad op = 1-09 using 
the methodology as described by lLoeb fc Gaudil (120031 ) . For Keplerian circular orbits, the 
(non-relativistic) radial velocity semi-amplitude is defined as 



K 



2nG 



1/3 



Mpsini 
r 2/3 



p J m: 

where G is the gravitational constant. The contribution of flux from Doppler beaming is 

-^dop 



(6) 



v4 dop sin(27T0) 



(7) 



tot 



The flux variations due to reflected/emitted light from the companion can be approxi- 
mated by 



ref 



ag ( — ) m 



F to t \ a 

where Aq is the wavelength-dependent geometric albedo, and ^ is the phase function for a 
diffusely scattering Lambertian Sphere 



sine/) + (tt — (bjeose 



7Y 



(9) 
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Values of the model parame ters are derived using a Levenberg-Marquardt least-squares 



X 2 approach (IPress et al.l Il992l ). We fit for the orbital period (P), epoch (T^), impact 
parameter (b), scaled planet radius (R p /R*), geometric albedo (A G ), secondary eclipse depth 
(ED), radial velocity semi-amplitude (K), ellipsoidal variations (A;ii)> an d mean stellar 
density (p*). We then use a Markov Chain Monte Carlo (MCMC) method (eg. Ford 2005) 
to estimate values and uncertainties using this initial model solution to seed the runs. Four 
chains of 10 6 samples each are run and the first 25% are discarded to account for burn in, 
allowing the Markov chains to stabilize. The median values of the best-fit parameters for 
Kepler-41b are given in Table 1, along with the la (68.3%) confidence intervals. We find M p 
= 0.598±g;||| Mj, R p = 0.996lg;^ Rj, yielding a mean planet density of p p = 0.74+°;^ gm 
cm -3 . The best-fit amplitudes of the occultation depth, ellipsoidal variations, and variations 
in reflected/emitted light were found to be 60±9, 4.5+lg, and 37.4l| 6 ppm, respectively. 



4. Confirmation Method 

Our confirmation method involves two main steps: (1) We first simulate third light 
contamination scenarios by injecting a wide range of dilution factors into the full phase 
photometric model and re-fit each diluted light curve model to the photometry. The results 
from these model fits set limits on the stellar parameters of possible blends; (2) We next 
compare these stellar parameters with stellar evolution models to eliminate star/companion 
systems that are unphysical or inconsistent with a stellar blend. The goal is to determine 
the probability that the planet-like signature could be caused by a contaminating star in the 
aperture. 

To model a full range of stellar blends, we iteratively dilute the best-fit model light 
curve with a dilution factor D which ranges from 1% - 100% of the transit depth (using 
1% intervals). For each value of D, we fit the diluted model to the light curve data and 
recompute x 2 (the goodness-of-fit estimator), p*, P, b, Rp/R*, ED, A e n, and A G . Figure 
4 shows results from these dilution fits for Kepler-41. The x 2 value and the change in x 2 
(Ax 2 ) as a function of injected dilution are shown in the top two panels. The lower six 
panels show the best fit results from six of the above parameters as a function of dilution. 
To determine the maximum third light from a potential blend, we solve for the dilution value 
at which Ax 2 increases by 1, 4, or 9, equivalent to la, 2a, or 3a (68.3%, 95.4%, or 99.7%) 
confidence intervals, respectively. These are shown by red, blue, and green horizontal lines 
in the top right panel of Figure 4. This resulted in maximum dilution values of 0.5 (la), 
0.6 (2a), and 0.67 (3c) (shown by the vertical red, blue, and green lines in the lower six 
panels of Figure 4). These constraints on the dilution values in turn place limits on the valid 
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parameter values of the companion which we can then compare to stellar evolution models 
to begin ruling out stellar blends. 



We use Yonsei-Yale (YY) stellar evolution models ( lYi et al.ll2003l ; iDemarque et al.ll2Q04i ) 
which provide a stellar age, T eSt , R*, p*, and logg for a given M* and Z. We start with a 
grid of stellar masses (M* = 0.4 - 5 M Q with 0.1 M Q increments) and metallicities [Z = 
0.00001, 0.0001, 0.0004, 0.001, 0.004, 0.007, 0.01, 0.02, 0.04, 0.06, 0.08) that covers the full 
range of input values to the YY models. For each (M m , Z) pair, we extract all YY models 
that have p* values within the constraints set by the dilution fits. For Kepler-41, the valid 
values of p* range from 1.17 gm cm -3 to a maximum value of 1.7795, 1.7800, or 1.7863 gm 
cm -3 for the la, 2a, and 3a constraints, respectively. Note the confidence intervals that 
provide constraints on the fit parameters are measured with respect to A% 2 rather than the 
distribution of each fit parameter. We require that the age from each model is less than the 
age of the Universe (taken to be 14 Gyr), and consider only those with ct-enhanced mixture 
equal to Solar mixture. The total number of models extracted for Kepler-41 was 1931, 1947, 
and 1950 for the la, 2a, and 3a cases, respectively. 

For each model, we derive the planet mass M p , radius i? p , equilibrium temperature T eqp , 
effective temperature T e fj p , and mean density p p as follows. We first take the p* value for 
the given model and interpolate over the valid range of p* (as shown in Figure 4) to find the 
corresponding dilution factor. This dilution value is then used to determine the values of the 
additional fit parameters (R p /R*, A e n, P, b, and ED) for that model. From these estimates, 
we can derive the remaining parameters that we use to characterize the star-companion 
system. 

The planet radius is found by R p = {R v l R*)R*. The planet mass is roughly proportional 
to the amplitude of the ellipsoidal variations (jPfahl et al.l 120081 ) and can be estimated using 
Equation 2. We solve for a/R* by combining the mean stellar density p* = M ;f /((4/3)7r_R 3 ) 
with Kepler's third law (for M p < M») M* = 4n 2 a 3 /(P 2 G), yielding 



a 

R* 



P*GP r < 
3vr 



(10) 



The inclination is computed from the impact parameter, which (for a circular orbit) is 
b = a/R*cosi. The planet mass is then 

(a/R*Y A e uM* 



(11) 
a eU sin i 

which, along with R p , provides a measurement of the planet density pp. We note that the 
amplitude of the Doppler beaming, if d etected (which is not the case for Kepler-4 1), can also 
be used to estimate the planet's mass (jShporer et al.l 1201 ll ; Barclay et al.ll2012bl ). 
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The planet equilibrium temperature can be estimated by 

(12) 

where A-q is the wavelength-integrated Bond Albedo (we use = 0.02, which is the ap- 
proximate value for hot Jupiters in t he Kepler bandpas s), and / is a circularization factor 




which equals 1 for isotropic emission ORowe et al.ll2008al ). 



The secondary eclipse depth is approximately equal to the ratio of the planet and star 
luminosities, ED fh L p /L*, and can be used to estimate the planet effective temperature 
T e g . Assuming blackbody radiation, the bolometric luminosity (the total amount of energy 
emitted across all wavelengths) of an object can be computed from the Stefan-Boltzmann 
equation L oc -R 2 T e 4 ff , where R and T eS are the radius and effective temperature of the body, 
respectively. The planet effective temperature can be solved from 

Rp\ 2 { T eKp 



ED - [t) \£) (13) 

Alternatively, we can compute L p /L* by integrating the planet and star Planck functions 
over the Kepler bandpass (A ~ 400 - 900 nm). The wavelength-dependent luminosity can 
be solved from 



L\d\ = Air R B\d\ (14) 



where B\ is the Planck function: 



B - 2hc2/X " (15) 

Here, h is Planck's constant, k is Boltzmann's constant and c is the speed of light. This 
latter method provides a more accurate estimate of T effp , and results using both methods to 
compute T e ff will be discussed in the next section. 



4.1. Results 

Our dilution models at the la, 2a, and 3a confidence levels combined with stellar 
evolution models yield a total of 5828 star/companion model configurations, each providing 
an estimate of stellar age, T cSf , R*, p*, log p, M p , _R P , T effp , T eqp , and p p . The next step is to 
examine these models and eliminate those that furnish unphysical stellar properties. 
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Figure 5 shows function of T e ff t for all available YY evolution models (shown by 

the red curves in each panel). Also shown are the values of mean planet density as a function 
of planet effective temperature (p p versus T eSp ) for each 3a dilution model (shown by the 
multi-colored tracks in each panel of Figure 5, each color representing a value of metallicity). 
In order to comprise a viable stellar blend, a dilution model needs to reside in a region that 
overlaps with a stellar evolution track. In the left panel of Figure 5, T c g p was calculated 
by integrating the Planck function over the Kepler bandpass (as described in the previous 
section). None of the resulting dilution models lie in the vicinity of any stellar evolution 
track, thus eliminating all potential stellar blends in the Kepler-41 photometry. Using this 
method to calculate T e g p , we can conclude that the companion to Kepler-41b is a planet. 

The right panel of Figure 5 shows the dilution models computed using values of T eSp that 
were estimated from the bolometric luminosities of the star and companion (which are not 
as precise as those computed by integrating the star and companion Planck functions, but 
are simpler to derive). In this case, a subset of the models do overlap with stellar evolution 
tracks. Although we have shown that we can use the alternative method of computing T e s 
to confirm the planetary nature of Kepler-41b, we show here how we can further examine 
these overlapping models to exclude them as potential third light contaminants (which may 
be necessary for confirming other planet candidates). 

Figure 6 shows T eqp as a function of T e g p for the subset of dilution models that are 
consistent with stellar evolution tracks (those shown in the right panel of Figure 5). For 
the companion to be of planetary nature, we expect a near balance in these temperatures, 
T e ff p ~ Tgqp, meaning that any incident energy upon the planet is re-radiated. If T e g p » 
T eqp , however, the companion must be burning Hydrogen, i.e., of stellar nature (although 
it is feasible that the object could be a young planet). For the case of Kepler-41b, T effp for 
these remaining dilution models aren't substantially greater than their corresponding values 
of T eqp . We therefore can't definitively rule out stellar blends using this comparison, but it 
may be useful for other planet candidate systems. 

We next compared the companion mass to the stellar mass (Figure 7) for the same 
dilution models as shown in Figure 6. The values of M* are all between 0.8 -1.1 M & 
whereas the dilution models all have masses below 0.004 M , well below the ~0.08 M mass 



limit required for Hydrogen burning (IKumari Il963l ). i.e., the companion cannot be a star. 



With this comparison, we can eliminate the remaining dilution models since we have shown 
that we cannot produce a proper stellar blend of any kind. 

Our future plans include using this new method of combining phase curve modeling 
with stellar evolution models to both confirm and characterize additional Kepler planets. 
The potential to detect occultations in planet candidate light curves can be determined by 
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combining signal-to-noise measurements with an assumption of an albedo of approximately 
30% (Rowe et al. 2013). Based on the planet and star characteristics tables from Batalha 
et al. (2012), we expect about two dozen planet candidates in the Kepler Field of View will 
have the potential to be confirmed with this method. 



5. The Albedo of Kepler-41b 



5.1. Hot Jupiters Albedos 

To date, Kepler's unprecedented photometric precision has enabled th e detection of 
planetary occultations in visible wave lengths for eight hot Jupiters: HAT-P-7b ([Christiansen et al 



2010 



Welsh et al. 



2010b 



), Kepler-5b 



Kipping fc Bakodl201ll;lDesert et al. 



2011a). Kepler-6b 



(IKipping &: Bakos 



2011 



Desert et al 



2011at). TrES-2b (IKipping fc Spiegel 



(IDesert et al.ll2011bl ; iBonomo et al.ll2012l ) and KOI-196b flSanterne et aLll2011bl ). 



2011af) . Kepler- 7b (IKipping fc Bakos 



2011 



Demory et al. 



201l|), Keoler-12b JFortnev et al.l boilh Kepler-17b 



The Spitzer Space Telescope gathered thermal planetary emission measurements in in- 
frared wavelengths for several dozens of hot Jupiters. The aforementioned Kepler detections 
allow us to probe giant planet irradiated atmosphere properties at optical depths that were 
not e xplored before, thereby constr aining further the energy budget of hot Jupiter planets 



[e.g., 



Madhusudhan fc Seagerl 120091 ). 



In our solar system, gas giant geometri c albedos range fr om 0.32 for Uranus to 0.50 for 
Jupiter in a bandpass similar to Kepler's (IKarkoschkal Il994l ) . This is mainly due to their 
low equilibrium temperatures, as compared to hot Jupiters, which allow the formation of 
cloud decks made of a mmonia and water ic e in their atmosphere that are highly reflective 
in visible wavelengths ( iDemory et al.ll2011bl ). 



Hot Jupiters emit very little in visible wavelengths. The albedos of hot Jupiters were 
expected to be low due to e f ficient reprocessing of stellar incident radiatio n into thermal 
emission ( IMarley et al.l Il999l ; ISeager fc Sasselovl |2000| ; ISudarsky et al.l 120031 ). In addition, 
the presence of alkali metals in hot Jupiter atmospheres (Na and K) as well as TiO and 
VO (at the hotter range) is expected to cause significant absorption at visible wavelengths, 
rendering most hot Jupiters dark. 

The first const raint on a hot Jupiter vi sible flux was obta ined with the MOST satellite 
(jWalker et al.l 120031 ) observing HD209458b (IRowe et al.ll2008bl ). The corresponding geomet- 
ric albedo 3a upper-limit of A g < 0.08 confirmed these earlier theoretical predictions. The 
majority of hot Jupiter occultations measured by Kepler photometry corroborate today the 
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hypothesis that hot Jupiters emit very little in visible wavelengths, their measured geomet- 
ric albedo being attributed to thermal emission leaking into shorter wavelengths rather than 
contribution from Rayleigh scattering, clouds or hazes. 

Remarkably, a few irradiated giant planets exhibit visible flux in the Kepler bandpass 
that exceeds the expected contribution from thermal emission alone. A recent detailed 
analysis of Kepler- 7b occultation measurements showed a significant departure of the mea- 
sured brightness temperature as compared to the equilibrium temper ature, suggesting tha t 
the planetary flux is dominated by Rayleigh scattering and/or hazes ( IDemory et al.ll2011al ). 
In addition, combining visible and Spitzer infrared occultation measurements showed that 
Kepler-12b also exhibits an excess of flux in the visib le, possibly indicatin g a reflective com- 
ponent in this low-density hot Jupiter atmosphere (IFortney et al.l 120111 ). Ideas that have 
been invoked to explain the wide variation in obse rved hot Jupiter albedos include vari- 
ations in planetary densities (Sudars kv et al.l 120031) and condensates phase transitions at 



narrow temperature ranges (IDemory et al.l 



2011a 



Kane fc Gelinol 120121 1 . 



5.2. Kepler-41b as another outlier? 

Our global analysis yields an occultation depth of 60 ± 9 ppm, which translates to a 
geometric albedo of A g = 0.23 ± 0.05. Using a blackbody spectrum for the host star, the 
corresponding brightness temperature is 2420 K, which is ~400 K larger than the maximum 
planetary equilibrium temperature, assuming zero-Bond albedo and no stellar incident energy 
recirculation from the day hemisphere to the night hemisphere. Kepler-41b shows similar 
brightness temperature excess as Kepler- 7b, possibly suggesting contribution from Rayleigh 
scattering and/or hazes. 

The planetary phase modulation, caused by the combination between reflected light 
and thermal emission, has an amplitude that is ~la smaller than the occultation depth 
and slightly offset from the mid-occultation timing. At the high atmospheric pressures 
probed by Kepler (P~l bar), we would expect the even temperature across hemispheres to 
yield a phase curve exhibiting only nominal modulation. This result suggests that either 
atmosphere dynamics at depth are significantly more complex than this description or that 
the phase curve modulation is dominated by reflected light instead of thermal emission. 
Detailed modeling and Spitzer infrared observations would be especially useful toward a 
precise constraint on the planetary energy budget and could unambiguously disentangle the 
thermal emission and reflected light components. 
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6. Summary 

We have presented a new method to confirm giant planets purely by analysis of the 
photometric light curve combined with stellar evolution models. We have developed a full 
phase photometric model that includes both primary and secondary transits along with 
flux contributions from ellipsoidal variations, Doppler beaming, and reflected/emitted light. 
We inject a full range of dilution values into the model light curve to simulate third light 
contamination from stellar blends, and iteratively fit each diluted model light curve to the 
photometry. We then compare these fit results to stellar evolution models to determine if 
any set of diluted model parameters are valid (meaning the star and companion have masses, 
sizes and orbits that are consistent with a stellar evolution model) and match the shape of 
the photometric light curve. 

We applied this m ethod to Kepler-41, a G6V star with a recently confirmed giant planet 



(ISanterne et al.l 12011a ). using Kepler photometry taken during quarters Ql - Q8. The 
phased light curve shows a clear secondary occultation with a depth of 60±9 ppm. The 
phase of this occultation is near = 0.5, indicating the orbit of Kepler-41b is likely nearly 
circular. We detected flux variations due to reflected/emitted light from the planet (with 
an amplitude of 37.4j|; 6 ppm) and ellipsoidal variations (4.5jJ;g ppm), the latter of which 
enables us to estimate the mass and density of the planet. We did not detect variations 
due to Doppler beaming, but these measurements - if detected in the light curves of other 
planetary candidates - can also be used to measure the planet mass. 

To determine whether any dilution models have properties consistent with a star, we 
first compared the denities and effective temperatures derived from the diluted models (p p 
versus T e fj p ) to all Yonsei-Yale evolution tracks (which provide p* as a function of T e g t for all 
valid stellar evolution models). To estimate T e s p for the dilution models, we computed the 
ratio of the planet luminosity to that of the star (L p /L*), which is approximately equal to 
the measured secondary eclipse depth, and solved for T eSp . We first computed L p /L* using 
the bolometric luminosities (the total amount of energy emitted across all wavelengths) for 
the planet and star. For comparison, we also computed these luminosities by integrating the 
Planck functions over the Kepler bandpass. This latter method to compute T e fj p resulted 
in unphysical star/companion parameters for all dilution models, thereby eliminating the 
possibility that the companion to Kepler-41b could be a stellar blend. Using values of T e g p 
that were computed from bolometric luminosities yielded a small subset of dilution models 
that were consistent with stellar evolution tracks. For these models, we further examined the 
temperatures and masses of each system to filter out additional inconsistencies. We found 
that all companion masses from these remaining diluted models were well below the ~0.08 
M & limit for Hydrogen-burning, indicating that the companion cannot be a star. Although 
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both methods to compute T e g p provided enough information to rule out stellar blends in the 
Kepler-41 photometry, we recommend computing T c g with the more accurate method of 
integrating the star and planet luminosities over the Kepler bandpass. 

Our best-fit model of Kepler-41b yields M p = 0.598±g;||| Mj and R p = 0.996^;^ Rj. 
From our analysis of the phase curve combined with stellar evolution models we can therefore 
independently confirm that Kepler-41b is indeed a planet. This confirmation method can be 
applied to additional Kepler planet candidates that show a clear occultation in their light 
curve. 

This paper includes data collected by the Kepler mission. Funding for the Kepler mission 
is provided by the NASA Science Mission directorate. Some/all of the data presented in this 
paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is 
operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. Support for MAST for non-HST data is provided by the NASA Office 
of Space Science via grant NNX09AF08G and by other grants and contracts. 
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Table 1. Star and Planet Parameters for Kepler-41 



Parameters 


Med 


Std 


+ 1<T 


-la 


Af. (M ) 


1.044 


0.062 


0.059 


-0.068 


R* (Rq) 


1.065 


0.040 


0.042 


-0.042 


Z 


0.036 


0.011 


0.010 


-0.012 


Tea, (K) 


5617 


143 


126 


-177 


log L, (I,© ) 


0.007 


0.064 


0.062 


-0.069 


logs (cgs) 


4.398 


0.028 


0.027 


-0.032 


p~, (gem -3 ) 


1.213 


0.099 


0.081 


-0.113 


Period (days) 


1.8555577 


0.0000003 


0.0000003 


-0.0000004 


T (BJD-2454900) 


70.180267 


0.000067 


0.000065 


-0.000070 


b 


0.475 


0.048 


0.045 


-0.044 


Rp/ R* 


0.09601 


0.00070 


0.00070 


-0.00070 


a/K* 


6.04 


0.16 


0.14 


-0.18 


R P (Rj) 


0.9956 


0.0375 


0.0393 


-0.0401 


a (AU) 


0.02999 


0.00060 


0.00052 


-0.00071 


M p (Mj) 


0.5978 


0.7803 


0.3839 


-0.5978 


Pp (g/cm 3 ) 


0.74 


0.98 


0.48 


-0.74 


i (deg) 


85.50 


0.47 


0.40 


-0.50 


Duration (h) 


2.335 


0.015 


0.015 


-0.016 


Depth (ppm) 


10653 


13 


13 


-13 


Occultation depth (ppm) 


60 


11 


6 9 


-9 


Ellipsodial (ppm) 


4.5 


3.3 


2.8 


-3.8 


Reflected/Emitted (ppm) 


37.4 


6.3 


6.1 


-6.6 


K (m/s) 


94 


123 


61 


-94 
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Fig. 1. — Kepler Ql - Q8 PA photometry for Kepler-41 is shown in the top panel. The 
cotrended light curve (middle panel) shows mostly spot activity. The light curve was Fourier- 
filtered to remove stellar variability (lower panel). 
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Orbital Phase 

Fig. 2. — The detrended data for Kepler-41 phased to the orbital period and binned to 0.005 
in phase. The red-lined data are centered on the transit and show the full phase and the 
green-lined data are centered on the occupation and magnified (see top and right axes). The 
green curve fits for ellipsoidal variations, Doppler boosting, and reflected/emitted light. 
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Fig. 3. — The best-fit model for Kepler-41b phased to the orbital period and magnified to 
show the occultation. Our full phase photometric model includes flux variations induced 
by the companion that can be decomposed. These include Doppler beaming (blue dotted 
curve), ellipsoidal variations (green dashed curve) and reflected/emitted light (orange dot- 
dashed curve). The sum of these three effects is shown in red. Note we did not detect 
Doppler beaming in the light curve of Kepler-41, but we include a description of this effect 
in this article because it may be applicable to other planet candidates. 



- 19 - 







Fig. 4. — Results from our dilution model fits. The goodness-of-fit estimator x 2 is shown 
in the top left panel as a function of dilution values that were injected into the light curve. 
We solve for the maximum allowed dilution (i.e., the maximum amount of 3rd light from a 
potential blend) by measuring where A% 2 changes by 1, 4, or 9 (corresponding to la, 2a, or 
3a), as shown in the top right panel by the red, blue and green horizontal lines, respectively. 
The lower six panels show six of the fit parameters as a function of dilution, and the red, blue 
and green vertical lines determine their range of valid values as constrained by the dilution 
fits. Comparison of each valid dilution model to stellar evolution models rules out massive, 
stellar objects, confirming the planetary nature of Kepler-41b. 
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Fig. 5. — The mean stellar density p* is shown here as a function of T e ^ t for all available 
Yonsei-Yale stellar evolution tracks (shown by the red curves in each panel). The companion 
p p and T e g p from the dilution model fits are overplotted for a range of metallicities Z (colored 
points in each panel). The dilution models in the left panel were computed using estimates of 
T e g p that were computed by integrating the planet and star Planck functions over the Kepler 
bandpass and comparing the ratio of the resulting luminosities to the secondary eclipse 
depth ED. All dilution models in this case are inconsistent with any stellar blend (there 
is no overlap with the stellar evolution tracks), and we can conclude that the companion 
to Kepler-41 is a planet. In the right panel, the dilution models were computed using T e g 
values that were calculated from the ratio of the planet and star bolometric luminosities 
(over all wavelengths). This was done to determine if this simpler method (albeit not as 
precise) to compute T e g p is sufficient to rule out potential blends. In this case, a subset 
of dilution models overlap with stellar evolution tracks and therefore need to be examined 
further (see Figures 6 and 7) in order to rule out stellar blends. 
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Fig. 6. — For the remaining valid dilution models (those that have parameters that overlap 
with stellar evolution tracks as shown in Figure 5), the equilibrium temperatures can be 
compared to the effective temperatures. To be of stellar nature, the values of T e g p for each 
model would need to be much greater than the corresponding values of T eqp (indicating that 
the companion is burning Hydrogen). In this case, the temperatures are comparable and 
cannot be used to definitively rule out stellar blends, but this comparison may be useful to 
confirm other planet candidates. 
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Fig. 7. — For the remaining valid dilution models (as shown in Figure 6), the relation between 
each companion mass M p and the corresponding stellar mass M* is shown here. All dilution 
models have a companion mass less than that needed for Hydrogen burning (~0.08 M©), 
indicating that the companion cannot be a star. With this comparison, we can eliminate 
these remaining dilution models and conclude that the companion to Kepler-41b is a planet. 
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